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ABSTRACT 


The development of a satellite infrared technique for estimating convective and stratiform 
rainfall and its application in studying the diurnal variability of rainfall in Amazonia are 
presented. The Convective-Stratiform Technique, calibrated by coincident, physically retrieved 
rain rates from the TRMM Microwave Imager (TMI), is applied during January to April 1999 
over northern South America. The diurnal cycle of rainfall, as well as the division between 
convective and stratiform rainfall is presented. Results compare well (a one-hour lag) with the 
diurnal cycle derived from TOGA radar-estimated rainfall in Rondonia. The satellite estimates 
reveal that the convective rain constitutes, in the mean, 24% of the rain area while accounting for 
67% of the rain volume. The effects of geography (rivers, lakes, coasts) and topography on the 
diurnal cycle of convection are examined. In particular, the Amazon River, downstream of 
Manaus, is shown to both enhance early morning rainfall and inhibit afternoon convection. 
Monthly estimates from this technique, dubbed CST/TMI, are verified over a dense rain gage 
network in the state of Ceara, in northeast Brazil. The CST/TMI showed a high bias equal to 
+33% of the gage mean, indicating that possibly the TMI estimates alone are also high. The root 
mean square difference (after removal of the bias) equaled 36.6% of the gage mean. The 


correlation coefficient was 0.77 based on 72 station-months. 


1 . 


Introduction 


In this paper, we present the application of and results from a satellite infrared (IR) technique 
for estimating rainfall over northern South America. Our objectives are to examine the diurnal 
variability of rainfall and to investigate the relative contributions from the convective and 
stratiform components. There is a paucity of hourly-reporting rainfall stations in northern South 
America in general and the Amazon Basin in particular. Hence, most of the literature describing 
the diurnal cycle of rainfall in this region has necessarily used satellite observations as a proxy 
for rainfall. One notable exception is Kousky (1980), who examined the diurnal cycle of rainfall 
in northeast Brazil from 10 years of gage observations. Estimates from IR-only methods include 
Meisner and Arkin ( 1987), who found a pronounced maximum in IR cloudiness at 1 800 LT over 
the interior of summertime South America. Janowiak and Arkin (1991) examined the period 
1986-1989 for seasonal and interannual variations in estimated rainfall. Horel et al. (1989) 
looked at 15 years of twice-daily out-going longwave radiation measurements to describe the 
annual cycle of convection in the Amazon Basin. More recently, Garreaud and Wallace (1997) 
examined 9 years of high spatial resolution IR data to define the diurnal march of deep 
convection over the tropical Americas and surrounding waters. 

With the launch of the Special Sensor Microwave Imager (SSM/1) in 1987, microwave (MW) 
algorithms were developed and applied. Examples include Negri et al. (1994), a 5-year 
climatology of warm-season tropical rainfall, and Negri et al. (2000), and a 10-year climatology 
over northern South America. Using four SSM/I observations per day, they showed dramatic 


differences in afternoon versus morning rain estimates. Methodology using a combination of 
MW and IR include Huffman et al. (1995, 1997) who made global, monthly rain estimates from 
a merger of IR, MW, gage and models. Other techniques adjust existing IR techniques by the 
MW measurements. Examples include Negri et al. (1993), a 3-year, warm-season study over 
Mexico, and Xu et al. (1999, 2000) who adjust the GOES Precipitation Index (GPI) by MW 
observations from SSM/I in an attempt to improve monthly rainfall estimates. Recently, the 
Passive Microwave and Infrared Algorithm (MIRA) [Todd et al., 2000] attempts to estimate 
small-scale rainfall. In the forerunner of the research described in this paper, Anagnostou et al. 
(1999) trained IR observations with three months of SSM/I observations to deduce the diurnal 
cycle of rainfall over northern South America. 

In this study, we apply the Convective-Stratiform Technique (CST) of Adler and Negri, 
[1988], The parameters of the original technique were recalibrated using coincident rainfall 
estimates derived from the Tropical Rain Measuring Mission (TRMM) Microwave Imager 
(TMI) and GOES IR ( 1 1 pm) observations. Henceforth, we refer to this new technique as the 
CST/TMI. 


2.0 Data and Methodology 

In this study, we apply the Convective-Stratiform Technique (CST) of Adler and Negri, 
[1988]. In short, the CST locates, in an array of IR data, all local minima in the brightness 
temperature field. After an empirical screening to remove thin, non-raining cirrus, the remaining 
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minima are assumed centers of convective rain. A stratiform rain assignment, based on the mode 
temperature of each cloud system, completes the rain estimation. The parameters of the original 
technique were recalibrated using coincident rainfall estimates derived from the application of 
the Goddard Profiling (GPROF) algorithm to TMI brightness temperatures [Kummerow and 
Giglio, 1994; Kummerow et al . , 1996]. We use an improved algorithm for convective/stratiform 
division in the GPROF as described by Olson et al., [2000]. 

Over northern South America, the GOES-east observes every 30 min; therefore, we were 
able to match a TMI image with an IR image within a 15 min window. Calibration was 
performed over land only, from 12N to 18S and from 82W to 34W. The TMI rain estimates 
contain the distribution of total surface rainfall as well as the classification of convective and 
stratiform rain. The matched data (165 TRMM overpasses) are in the period from January to 
April, 1999, the first two months of which coincide with the period of the Large-Scale 
Biosphere-Atmosphere (LB A) experiment in Brazil. 

The goal of the recalibration is to determine parameters so that the CST will be able to 
reproduce the total rain volume, total rain area, and the TMI-observed division between 
convective and stratiform rain. We adjust the CST on a statistical basis, over the 4-month 
calibration period, without explicit constraints on instantaneous rain estimates. The recalibration 
takes the following form: 
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a) determination of convective cores using a discriminant analysis. This is based on the 
identification of points of local IR temperature minimum (T mjn ) and the computation of the 
deviation of the T min from the background temperature (defined as the average temperature of the 
eight surrounding pixels). Figure 1 shows the analysis of T min in the 2-dimensional space defined 
by IR temperature and the deviation of the T min from the background. Points are defined by the 
corresponding TMI-retrieved rain rate as either convective (C) or non-convective (N). Plotted 
data represents a 1% random subset of the total dataset used to derive the discriminant line. T min 
points above and to the left of the discriminant lines are deemed convective. The resultant 

discriminant lines are defined by: 

1.25- T m , n - 3.16-Deviation < 254.7 and Deviation < 2.23 

b) determination of the convective rain area assigned to each T min . This is based on the 
computation of a, the constant of proportionality between the total convective area in the TMI 
estimates and the parameter X(T cU)uJ - T mm ). For a given convective T min , a • (T cloud - T min ) is the 
number of pixels to which we assign a convective rain rate centered on the T min . In this study 
T C ioud is a constant set to 253 K, a commonly used threshold for convective clouds. Based on the 
rain area determined from the TRMM convective area estimates, a = 0.61. For T min = 203 for 
example, 32 pixels (4 by 4 km) would be assigned the convective rain rate. 

c) determination of the mean convective rain rate. We set the convective rain rate equal 
to the conditional rain rate of the TMI/GPROF estimates, namely 18.9 mm h . This allows the 
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total convective rain volume derived from CST /TMI to equal that calculated from the TMI 


alone. 

d) determination of the threshold temperature to define the stratiform rain area. We 
choose a threshold such that the total cloud area colder than this threshold (excluding those 
pixels already assigned convective rain) is equal to the stratiform rain area derived by 
TMI/GPROF. This threshold, calculating using cumulative area matching, was found to be 219 
K. 

e) determination of the mean stratiform rain rate. As in the case of determining the 
convective mean rain rate, the conditional TMI/GPROF stratiform rain rate of 2.6 mm h was 
used. This forces the total stratiform rain volume from CST/TMI to equal that of the TMI alone, 
integrated over the 4-month calibration period. 

In summary, the CST/TMI finds local IR minima, decides if they are convective features, 
then assigns a rain area and rainrate based on the TMI calibration parameters. A stratiform rain 
area is then defined, and a lower rain rate is assigned to those cloud pixels colder than the 
stratiform IR threshold and not previously assigned convective rain. We have chosen two rain 
rates to make the technique simple and efficient. An example of the instantaneous application of 
this technique is shown in Figure 2. The four panels indicate (clockwise from the upper left), a 
GOES IR image of southwestern Brazil; the CST/TMI rainfall estimate; the rain rate from the 
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TRMM radar algorithm 2A25; and the GPROF-retrieved rain rate applied to the TMI brightness 
temperatures. For reference, the outline of the TMI swath appears on all images. 

3. Results 

3.1 The Diurnal Cycle in the TOGA Radar Area 

Rain estimates over northern South America (12N to 18S, 82W to 34W) were made from 30 
min interval GOES data for the period January-April, 1999. The first two months of this period 
coincided with the LBA experiment in Rondonia, (southwestern Brazil). TOGA radar data 
provided a comparison between the mean radar rain rates and convective/stratiform division with 
those estimated by the CST. The satellite estimates were matched in both time (5 min) and space 
(accounting for bad and missing data) to the gage-calibrated radar estimates [Anagnoustou et al., 
2000], Figure 3 (left) displays the mean rain rate (mm h ') from the TOGA radar at 1 km height 
for the period 7 Jan to 28 Feb 1999. The asterisks mark the location of the four main clusters of 
rain gages. Problems with beam blockage and ground clutter along several azimuths necessitated 
the creation of a quality control mask (Figure 3, right). This mask represents the frequency of 
occurrence of rainfall in the interval 7-13%. (Higher than 13% was determined to be ground 
clutter; lower than 7% was beam blockage). Instantaneous radar estimates were averaged only 
within the confines of this mask, and it is apparent that problems remain with the radar data. 

Figure 4 shows the diurnal variation of four rain parameters during the experiment period of 
7 Jan to 28 Feb 1999. The upper left panel of the figure displays the unconditional rain rate for 
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the CST/TMI estimates (green), TOGA radar (red) and rain gages (blue). The upper right panel 
shows the percent of the radar area covered by rain greater than 1 mm h 1 . The lower left panel 
tracks the variation of the convective rain rate (the percent of the rain volume due to convective 
rain). Finally, the lower right panel reveals the variation of the percent of the rain area due to 
convective rain. The satellite estimates of the mean rain rate (top left) are in good agreement 
with both the radar and gages. They lag the radar (and gage) estimates by one hour, essentially 
missing the earliest stages of convection. By comparison, results using a simple threshold 
method (the GOES Precipitation Index or GPI, not plotted) showed a three-hour lag in the 
maximum rainfall. The mean values of the satellite, gage and radar estimates are, respectively, 
0.35, 0.35 and 0.43 mm h' 1 . Compared to the radar estimates, the CST/TMI overestimates the 
fractional convective rain volume (bottom left) by a factor of 1.36, while underestimating the 
fractional convective area (bottom right) by a factor of 1.57. An underestimation of the total rain 
area by a factor of 1.52 is also noted. Despite the biases, the phase of all four of the satellite- 
estimated parameters seems to well represent (within an hour of the maximum), the phase of the 
radar-derived rainfall. The satellite estimates reveal that the convective rain constitutes, in the 
mean, 24% of the rain area while accounting for 67% of the rain volume. 

3.2 Results for Northern South America 

Figure 5 shows the application of the CST/TMI to 30 min interval data from January-April 
1999 over the northern portion of South America. Results (in mm h ') are binned into hourly 
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categories and adjusted for local time (3 zones). The final panel displays the mean rain rate for 
the entire 4-month period. Among the interesting features revealed in this presentation are: 

• The development of afternoon (12-20 LT) squall lines along the Northeast Brazilian 
coast. Such squall lines were studied by Garstang et al . , [1990]. 

• A morning maximum and afternoon minimum in rainfall along the Amazon River east of 
Manaus. This is described in more detail in the next section; 

• A morning rain maximum in the Gulf of Panama, no doubt in response to a land breeze 
along the highly convergent coastline; 

• The onset of deep convection in the Amazon Basin around 13 LT, maximizing at 15-16 
LT; 

• A late evening (23 LT) to early-morning (05 LT) maximum along the eastern slopes of 
the southern Andes Mountains, possibly the result of a mountain/valley circulation; 

• An afternoon (13-18 LT) maximum along the western slopes of the southern Andes. 

Many more small-scale features are revealed when these 24 images are viewed in time- 
sequence. A QuickTime'"' movie of this sequence is available by contacting the first author. An 
alternate method of displaying these hourly data is shown in Figure 6, the arithmetic difference 
of the CST/TMI estimates at 18 LT and 06 LT. These times correspond with the overpass times 
of the Special Sensor Microwave Imager (SSM/I). Areas with an afternoon maximum are 
displayed in yellow to red shades, while areas of morning maximum appear blue to purple. In an 


8 



effort to validate these results, we present Figure 7, the estimated rainfall difference (18 -06 LT) 
derived from the application of the GPROF algorithm to SSM/I brightness temperatures from the 
F13 DMSP satellite. Spatial resolution is 0.5° and the temporal sampling is about 20 
observations/month at both 18 and 06 LT. Note the similarity of features, particularly along the 
coasts and rivers, which lend credibility to the IR estimates. 

3.3 Effects of the Amazon River on the Rainfall 

Figure 8 is a blow-up of the mean hourly estimates in Figure 5, along the confluence of the 
Negro and Solimbes Rivers near Manaus, Amazonia, Brazil. Here we find a late-evening (00 
LT) rainfall maximum that elongates and moves westward along the river until about 10 LT. 
Convection begins by 12 LT, and there is a pronounced minimum in rainfall along the river from 
about 15 to 18 LT. 

Figure 9 is a time/latitude display of the mean estimated rain rate for January- April, 1999 
(left) and January-April, 2000 (right) for a cross-section through the Amazon River at 56W. The 
mean position of the Amazon River at this longitude is noted. Downstream of where the Negro 
and Solimbes Rivers merge, we find an early morning (3 LT) maximum along the river. Rainfall 
avoids the river in the afternoon (12 LT and later). The river seems to be generating a land/river 
circulation, which enhances early morning rainfall but inhibits afternoon rainfall along the river. 
Such phenomenon have been studied by Oliveira & Fitzjarrald, [1993]. The pattern is 
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duplicated a year later, January-April, 2000 (right), demonstrating the repeatability of the 
phenomenon. 

3.4 Diurnal curves for other regions 

Buoyed by the success of the technique in an area with ground-truth, we examine the diurnal 
variation of rainfall in other regions. In Figure 10, we display a time series of the estimated total, 
convective and stratiform rain rates (mm/h) for three regions of interest. In the region 
encompassing much of the Amazon Basin (top), peak convective rainfall occurs at 15-16 LT. 
The convective rainfall comprises about 74% of the total rainfall at the time of the maximum, 
and 67% overall. The stratiform estimates lag the convective by two hours. Estimates for a 
smaller subset of the Amazon Basin, along the Amazon R. (downstream of the confluence of the 
Negro and Solimoes Rivers) are shown in the middle panel of Fig. 3. A pronounced rainfall 
maximum at 3 LT is noted. The bottom panel of Fig. 3 shows the diurnal curve for the northeast 
Brazil area, where coastal squall lines form almost daily. The peak estimated rainfall is 1.3 
mm/h at 17 LT, and at the maximum, has a percentage of convective rainfall of 65%. The 
stratiform lags the convective by 2 h. 

3.5 Monthly Verification in the State of Ceara, Brazil 

Figure 1 1 is an example of the monthly rainfall (April 2000) in the state of Ceara in northeast 

Brazil. Monthly rainfall is highest along the coast, and tapers off as one moves inland, until 
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higher elevations are reached in the southern part of the state. This gage data was averaged in 
the nine 1° squares outlined in the figure for comparison to the area-averaged CST/TMI 
estimates. Figure 12 shows a scatterplot of the estimates versus the gages for the wet season for 
two years (January- April, 1999-2000). Standard deviations of the mean gage rainfall for the nine 
grid squares in Fig. 1 1 are plotted to demonstrate the high spatial variability of rainfall. The 
CST/TMI shows a bias (+33% of the gage mean). The root mean square difference after 
removal of the bias (RMSD-br) is 36.6% of the gage mean, but the technique is not specifically 
calibrated in this region. We are encouraged by the high correlation coefficient of 0.77. 

4. Conclusions 

In this paper, we have presented results from a satellite infrared technique for estimating 
rainfall over northern South America. The parameters of the technique were calibrated using 
coincident rainfall estimates derived from the application of the GPROF algorithm to TRMM 
Microwave Imager brightness temperatures. Our objectives were to examine the diurnal 
variability of rainfall at high space and time resolution, and to investigate the relative 
contributions from the convective and stratiform components. TOGA radar data during the LBA 
experiment provided radar ground truth for both the CST rainfall and convective/stratilorm 
estimates. 

The diurnal cycle of rainfall in the LBA radar area was constructed. The peak satellite- 
estimated rainfall lagged the radar by one hour. Identification of the onset of convection was the 
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primary cause of this lag. Estimates of the division of rainfall into convective and stratiform 
components were also accomplished. The phase of these components agreed with the radar 
estimates, though the convective volume (area) was overestimated (underestimated). The satellite 
estimates revealed that the convective rain constitutes, in the mean, 24% of the rain area while 
accounting for 67% of the rain volume. Local circulations were found to play a major role in 
modulating the rainfall and its diurnal cycle. These included land/sea circulations (notably along 
the northeast Brazilian coast and in the Gulf of Panama), mountain/valley circulations (along the 
Andes Mountains), and circulations associated with the presence of rivers. This last category 
was examined in detail along the Amazon R. east of Manaus. There we found an early morning 
rainfall maximum along the river (3 LT at 56W). Rainfall avoids the river in the atternoon ( 12 
LT and later. The width of the river seems to be generating a land/river circulation that enhances 
early morning rainfall but inhibits afternoon rainfall. 

Verification of the monthly CST/TMI results was performed using rain gages in the northeast 
Brazilian state of Ceara. The CST/TMI showed a high bias equal to +33% of the gage mean. 
The root mean square difference after removal of the bias was 36.6% of the gage mean. We were 
encouraged by the high correlation coefficient of 0.77 for a technique is not specifically 
calibrated in this region. 

Perhaps the most important future application of this technique may be assessing the impact 
of deforestation on the both the rainfall and it’s diurnal cycle. We intend to run the CST/TMI for 
a full year, and compare monthly and seasonal total and hourly rain estimates with maps of 
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forested and deforested regions of Amazonia. We hope that circulations like that described over 
the Amazon River will exist over these regions, and manifest themselves in the satellite rainfall 
estimates. A second application involves the application of this technique to global, high- 
resolution, geosynchronous IR data sets ( Janowiak et ai, 2001). 
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Figure Captions 


Figure 1. Discriminant analysis of IR temperature minima (T min ) in the 2-dimensional space 
defined by IR brightness temperature and the deviation of the T mjn from the background. 
Points are defined by the corresponding TMI-retrieved rain rate as either convective (C) or 
non-con vective (N) rain. Plotted data represents a 1% random subset of the total dataset used 
to derive the discriminant line. 

Figure 2. An example of the application of the CST/TMI to one IR image. The outline of the 
TMI swath is overlain on all images. 

(Upper left): GOES IR image of southwestern Brazil scanned at 2358 GMT 22 Jan 1999; 
(Upper right): CST/TMI rainfall estimate. Convective cores (green) have a mean rain rate of 
18.9 mm h 1 while stratiform rain areas (purple) have a mean rain rate of 2.6 mm h 
(Lower left): Rain rate retrieved by application of the GPROF algorithm to the TMI 
microwave brightness temperatures. Overpass occurred 7 min after the time of the GOES 
image; 

(Lower right): Rain rate from the TRMM radar algorithm 2A25. 

Figure 3. TOGA radar at 1 km height for the period 7 Jan to 28 Feb 1999. The asterisks mark 
the locations of the four main clusters of rain gages. 

(Left): Mean rain rate (mm h' 1 ) 
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(Right): Quality control mask applied to account for missing or bad azimuths and ground 


clutter. 

Figure 4. The diurnal variation of four rain parameters in the period 7 Jan to 28 Feb 1999. 

(Upper left): The unconditional rain rate for the CST/TMI estimates (green), TOGA radar 
(red) and rain gages (blue). 

(Upper right): The variation the fraction of the radar area covered by rain greater than 1 mm 
h 1 ). 

(Lower left): The variation of the percent of the rain volume due to convective rain. 

(Lower right): The variation of the percent of the rain area due to convective rain. 

Figure 5. Application of the CST/TMI to 30 min interval data from January-April 1999. Results 
(mm h ’) are binned into hourly categories and adjusted for local time (3 zones). The final 
panel displays the mean rain rate for the entire 4-month period. 

Figure 6. The arithmetic difference of the CST/TMI estimates at 18 LT and 06 LT for the 
period January-April, 1999. Areas with an afternoon maximum are displayed in yellow to red 
shades, while areas of morning maximum appear blue to purple. Units are mm h '. 

Figure 7. Same as Figure 4 except the rainfall estimates are derived from the application of the 
GPROF algorithm to SSM/I brightness temperatures from the F13 satellite. Spatial 
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resolution is 0.5° and the temporal sampling is about 20 observations/month at both 18 and 


06 LT. 


Figure 8. A blow-up of the mean hourly estimates (mm h ') in Figure 3, along the confluence of 
the Negro and Solimoes Rivers near Manaus, Amazonia, Brazil. 

Figure 9. A time/latitude display of the mean estimated rain rate (mm h 1 ) for a cross-section 
through the Amazon River at 56W. 

(Left): January- April, 1999 
(Right): January - April 2000 

Figure 10. The diurnal cycle of CST/TMI rain estimates for three regions of South America: 
(Top): the Amazon Basin (0-10 S, 50-75W) 

(Middle): the Amazon River (2-3 S, 55-59 W) 

(Bottom): the northeast Brazilian coast (2-5 S, 42-50 W) 

Figure 1 1. Monthly gage rainfall for April 2000 in the state of Ceara in northeast Brazil. 

Figure 12. Scatterplot of CST/TMI versus gage rain estimates of monthly rainfall, averaged to 1° 
squares, for the 8 months January-April, 1999-2000. The vertical lines with bars represent 


20 



one standard deviation of the mean gage-averaged rainfall for the nine outlined 1° squares in 
Figure 1 1 . 
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Diurnal cycle in the Amazon Basin area (0-1 OS, 50— 75W) 
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Figure 10 





Ceara gages Apr 2000 
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Figure 11 


Raingage measured rainfall (m 


1 — deqree overoqed monthly roinfoll over Cearo, Brazil 
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Num = 72 

Gage mean = 151.2 
Sat. mean = 201.6 
Bias = 50.4 ( 33.3%) 
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Correl. = 0.766 
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Figure 12 



